Nanomechanical effects in an Andreev quantum dot 
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We consider a quantum dot with mechanical degrees of freedom which is coupled to supercon- 
ducting electrodes. A Josephson current is generated by applying a phase difference. In the absence 
of coupling to vibrations, this setup was previously proposed as a detector of magnetic flux and we 
wish here to address the effect of the phonon coupling to this detection scheme. We compute the 
charge on the quantum dot and determine its dependence on the phase difference in the presence 
of phonon coupling and Coulomb interaction. This allows to identify regions in parameter space 
with the highest charge to phase sensitivity, which are relevant for flux detection. Further insight 
about the interplay of such couplings and subsequent entanglement properties between electron and 
phonon degrees of freedom are gained by computing the von Neumann entropy. 

PACS numbers: 74.45.+C, 73.21. La, 74.78.Na, 46.40.-f 



I. INTRODUCTION 

The Josephson effect is one of the most striking man- 
ifestation of phase coherence in macroscopic objects. A 
nondissipative current^ can flow through a junction be- 
tween two superconductors, provided that there is a 
phase difference between them. Early Josephson junc- 
tions consisted of an oxide layer or a normal metal sand- 
wiched between the superconducting leads^i but progress 
in nanofabrication techniques have allowed to imbed 
mesoscopic devices into the junction. "^"^ One of the most 
commonly studied of such devices is the quantum dot. 
Quantum dots typically represent a normal metal island 
with resonant levels and possibly charging effects. In the 
context of Josephson transport, it has been shown the- 
oretically that the charge on such quantum dots can be 
tuned either by applying a gate voltage to the dot or 
by varying the phase difference between the supercon- 
ductors^ This continuous tuning of parameters allows 
the dot charge to deviate from an integer number. Of 
importance in such a system is that the tuning param- 
eters can trigger a transition of the ground-state from 
a singlet (zero or double electron occupancy with oppo- 
site spins) to a doublet (single electron occupancy with 
spin up or spin down).— lii In Ref. [Tl|, it has been pro- 
posed that the sensitivity of the dependence of the charge 
with respect to the flux could in principle be exploited to 
measure rather precisely the magnetic field in the loop, 
in the same spirit as a superconducting quantum interfer- 
ence device. The measurement of the charge itself could 
possibly be performed using a single-electron transistor 
coupled electrostatically to the dot in the junction. 

At the same time, in nowadays experiments, one has 
the possibility either to taylor artificial quantum dots 
and to embed them in a circuit, or alternatively to use 
existing nano-objects for the same purpose. In particu- 
lar, carbon nanotubes contacted to metallic or supercon- 



ducting leads^"= have been shown to behave like quantum 
dots, with the advantage that they can be influenced by 
nearby metallic gates .<^ There are also attempts to place 
single molecules in the junction between two reservoirsi^ 
In such systems, the vibrational degrees of freedom may 
affect electron transport in two ways. First, there are al- 
ways vibrational degrees of freedom associated with the 
material surrounding the molecular quantum dot. Such 
phonons typically constitute a source of relaxation and 
decoherence mechanism for quantum transport J^^— Sec- 
ond, the quantum dot itself may have internal vibrational 
degrees of freedom, which are coupled to the charge of the 
quantum dot<^"— We focus on the latter mechanism in 
this work. A number of previous works have addressed 
this issue for nonequilibrium transport with normal metal 
contacts)^"— or for the supercurrent through a vibrating 
nano-ob j ects ^— 

With this paper we want to address the issue of the 
phase sensitivity of the charge in an Andreev quantum 
dot, taking into account the presence of electron-phonon 
interaction. The goal is to determine the impact of the 
phonon coupling on the measurement scheme. Starting 
from a microscopic Hamiltonian model, we will compute 
the equilibrium properties of the system for various pa- 
rameters, in the regime where the superconducting gap 
is much larger than all other relevant energies in the sys- 
tem. 



II. MODEL 

Two typical setups are depicted in Fig. [1] On the right 
hand side is a generic setup where the central island is 
coupled to the right and left leads and which has a (sin- 
gle) vibrational degree of freedom. On the other hand, 
the setup on the left side [see Fig. HJa)] of the figure 
represents a single- wall carbon nanotube (SWNT) which 
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FIG. 1: (a) SWNT suspended between two superconducting 
leads with phase difTerence ip. The charge of the SWNT Q 
is displaced by an external electrical field E. Two additional 
gates form the electron's resonances between them, (b) Model 
setup represents a quantum dot with oscillator degree of free- 
dom coupled to superconductors through tunnel junctions. 



is suspended between two superconducting leads. Ad- 
ditional gates placed above the nanotube allow to define 
precisely the extent of the quantum dot, and therefore al- 
low to modulate its energy levels. An overall gate voltage 
allows to apply an electric field to the whole structure. In 
this second setup, several vibrational modes are known 
to exist 4^1^ 

We discuss primarily the case of a short junction (its 
length L is much less than superconducting coherence 
length ^) where the normal island can be described as a 
zero-dimensional (OD) object. Due to the shortness of the 
junction, in practice this setup effectively describes real 
OD molecular quantum dots. Such setups were realized 
experimentally'^'"' and described theoreticallyj^s^ii^ 

The model of the Andreev quantum dot is described 
by a total Hamiltonian which includes the dot and its 
internal degrees of freedom, the leads, and the tunnel 
coupling between the latter two. 



H — Ho + Hs + Ht^. 



(1) 



The first term Ho is the quantum dot, which contains for 
simplicity a single level and a discrete phonon spectrum, 

Ho = J2f^^Mh+ 

i 

[eo-^A.(6, +6l)] J2 {f^.-l)+Un^n^ (2) 

with ficr = d\d^\ dj., d^ are electron creation and an- 
nihilation operators for the dot. The energy e^ is the 
dot level, which can be tuned by a gate voltage, and 
which is measured with respect to the Fermi energy of 
the leads. Note that this simple model can represent a 
more realistic multi-level dot, when one level only con- 
tributes significantly to the electronic transport because 
the spacing between the dot levels is large compared to 
the superconducting gap and the coupling to phonons. 
Each term Mlib\h^ denotes the phonon energy of the vi- 
bration mode i in the dot {b\ and 6^ are phonon creation 
and annihilation operators), Ai is the electron-phonon 



coupling in this mode; index i runs over all mechan- 
ical modes i ~ 1, 2, . . . ,iVfnodcs- The electron-phonon 
coupling mechanism is described by terms XiE^ where 
Xi = y^h/2A'Iifli{b^ + bj) are displacements in an exter- 
nal electrical field E. The charge of the dot is attracted 
by external gate voltage, which leads to a change in its 
position. The deformation leads to the changing of the 
ground-state energy and, therefore, of the charge of the 
dot. 

In the sum J2cri''^<^~^/'^) ^^^ constant 1/2 is subtracted 
to "symmetrize" the matrix elements of the Hamiltonian. 
U describes the Coulomb interaction. The lead Hamil- 
tonian describes two BCS superconductors [with a lead 
index ^ = L,R (left, right)]. 
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(3) 



with an energy dispersion in superconducting leads ^k — 
h^k'^ /2m — Ep and an absolute value of the gap A in the 
bulk of the superconductors. The electron hopping term 
between dots and leads reads. 
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(4) 



where T^.r = ^L.RO'ze^*"'^'^/^ and ti's are tunneling am- 
plitudes between superconductors and the dot. ip is a 
superconducting phase difference. 

Calculations of observables for this system in ther- 
mal equilibrium typically start from the calculation of 
the partition function Z = Tr{exp(— /3iJ)}, where /3 = 
l/kuT is the inverse temperature. The Josephson cur- 
rent is then proportional to the logarithmic derivative 
with respect to the phase difference tp, and the charge on 
the dot is the derivative of the free energy with respect 
to the level position. In previous works using functional 
integral approaches''^ it was noted that because the to- 
tal Hamiltonian is quadratic in the lead fermion opera- 
tors, a partial trace over such degrees of freedom could 
be performed. This gives rise to an effective action with 
a dot fermion self-energy which contains retardation ef- 
fects, and which couples fermion operators of the same 
nature, but with opposite spins. This coupling is a man- 
ifestation of electron pairing phenomena at the level of 
the dot due to the proximity with the superconducting 
leads. In Ref. ISEl the calculation of the partial trace 
over the leads of the partition function was performed 
in a similar manner, nevertheless using an operator ap- 
proach. Furthermore, the assumption that \s-o\, U, F, 
hfl ^ A (the so-called A — > cx) limit) allowed there to 
neglect the retardation effect and to therefore derive an 
effective Hamiltonian for the dot-phonon system, 

1 
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+ r[did^ + h.c] + Uh^fii, (5) 



where f = Tcos{ip/2). The escape rate T = 27ri/(0)|ip 
(or resonance width of the dot), assuming a constant den- 
sity of states z^(0) near the Fermi energy of the metal 
in the normal state. We assume a symmetric setup 
I^Lp = I^Rp = |ip for the remainder of this study. The 
effective Hamiltonian Eq. ([5]) of the large A limit consti- 
tutes the starting point of our calculation. 

The matrix elements of the Hamiltonian with respect 
to the dot electron states \iy)ci are now computed. These 
states are: |0)ei (zero occupation), |t)ci = rf||0)oi, |i)ci 
= dt|0)oi (single occupation), |2)ei = (il(it|0)oi (double 
occupation), which means that from the electron point 
of view, the only off-diagonal part of H originates from 
the coupling to the leads and involve either zero or double 
occupancy states. 



H,,, = Y, f^Mh + k - 5] Kik + bl) 



diag{-l,0,0,l} + f 
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L/diag{0, 0,0,1}. 



(6) 



7^-1/2, 



and 



The basis of phonon states is |n)ph,i = (&])"(; 

|0)ph,i- Using relations i,ph("i|fo|&j|n)ph,i = n 

i,ph(w|6j -I- 6j|n)ph,i = \/n5m,n-i + V" + ^Sm,n+i, we 
can thus generate the phonon matrix elements of Eq. ^. 
For one single phonon mode the matrix representing the 
full Hamiltonian reads. 



H„ 



Gi Ai 

Ai 62 A2 

Ai 63 A3 

A3 64 



(7) 



We use greek indices for matrix elements in electron sub- 
space and latin ones for phonons. The 4x4 matrix blocks 
Qn and A„ are defined by 
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(8) 



and A^iy^„ = diag{— 1, 0, 0, Ij^/nA; they describe electron 
degrees of freedom with n phonons and electron-phonon 
coupling, respectively. The generalization to an arbitrary 
number of phonon modes can easily be obtained by mul- 
tiplication of Hilbert subspaces for each phonon modes. 
The eigenstates of the effective Hamiltonian can be cal- 
culated numerically by truncation of the matrix (trunca- 
tion of the number of phonon states). In practice, we 
took about 20 phonon states for A/0 = 3 and about 70 
states for A/0 = 5; these numbers are nearly independent 
of the number of modes A^modcs- 



III. STATES WITHOUT PHONONS: A = 

The electron states case were described using electron 
representation's and electron-hole Bogoliubov superposi- 
tion.— The effective Hamiltonian reduces to 4 x 4 matrix 
60 [sec Eq. (jH])] with singly degenerated (singlets) eigen- 
states 



E< 



0,2 



C//2T\/(^ 



C//2)2 + f 2 



and doubly degenerated (doublet) eigenstate 

El = E^^i = 0. 



(9) 



(10) 



In the absence of Coulomb interaction U = the eigen- 
values are ordered as Eq < Ei < E2 and the ground-state 
is always formed by the state with energy Eq (pure hole- 
like state in Bogoliubov representation). For U > the 
ground-state can be formed by the singlet |0)ei [the sin- 
glet region in 2D plane (eo, v)] or by the doublet |l)ei 
(doublet region), but never by |2)oi. For arbitrary finite 
U the doublet region exists if 

(eo + C//2)2 + f2< ([//2)2, (11) 

where F = T cos{(p/2). 



IV. REGIONS OF THE SINGLET AND 
DOUBLET STATES 

The existence of the doublet state as the ground-state 
of the system is important for this system and we dwell 
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FIG. 2: Singlet-doublet jump in the (<^, Ed) plane; the sin- 
glet regions lie outside the color loops, the doublet regions 
lies inside the loops, (a) For zero electron-phonon interaction 
A = and different Coulomb interaction fZ/F = (empty 
region), 2/3 (magenta; smallest region), 4/3 (blue), 2 (red; 
largest region). It is described by Eq. (Illf) . (b) For differ- 
ent phonon-electron interaction strength A = (red; largest 
region), A = 2 (cyan), A = 4 (green; smallest region) and 
fixed Coulomb interaction t//F = 2. The single-mode oscil- 
lator has a frequency which is much smaller than the tunnel 
resonance width Ml/T — 0.05. We see the interplay between 
Coulomb interaction and electron-phonon coupling; the first 
one expands the doublet region in the ip direction, the second 
one squeezes it in the Sd direction. 



on this more. For zero electron-phonon interaction A = 
the doublet region is specified by Eq. ([TT|) and its form 
is represented in Fig.[2lja). The highest value of U corre- 
sponds to the largest doublet region; with decreasing of 
U this region becomes smaller and smaller. At U = its 
disappears. 

If we start with some fixed finite U then the "area" of 
the doublet region decreases with electron-phonon cou- 
pling A. The evolution of the doublet region is plotted 
in Fig. EJb) for different values of the electron-phonon 
coupling constant. The largest loop corresponds to the 
smallest (zero) A. Upon switching A, the reduction of this 
region is barely noticeable, it acts mostly on the level po- 
sition range as it still approaches the phase values and 
2tt. There is a competition between charge repulsion ef- 
fects on the dot and the presence of the electron-phonon 
coupling, which can be understood to be playing the role 
of an effective attractive, retarded, interaction. This ex- 
plains the reduction of the doublet region. 

At horizontal line ip = tt the decreasing of the doublet 
region in e^ direction can be described by inequation 
\eo + U/2\ < U/2 — \^ /Ml. The nonzero electron-phonon 
coupling acts as the negative Coulomb interaction (in the 
sense of size of the doublet region) . It implies that for all 
A's larger than 



f/ = 



Ac = y/mu/2 



(12) 



the doublet region does not exist for any values of e^ and 

V. CHARGE OF THE ANDREEV DOT 

The charge Q of the nanotube/quantum dot [in a given 
quantum mechanical state, e.g., some eigenstate of ([5])] 
can be calculated by taking the derivative of its energy 
E (in this particular state) with respect to the external 
gate potential Vg (or Ed/g), 



Q = edE/d£j,. 



(13) 



If one is interested in the charge of the ground-state, then 
ground-state energy should be taken. 

The same result can be obtained by averaging the 
charge operator 



Q 






and the corresponding matrix elements 
Q^, = ediag{-1, 0,0,1} 



(14) 



(15) 



over the needed state [Eq. ((T5]) is written in the electron 
subspace; it should be multiplied by the unity matrix in 
the phonons subspace]. In what follows we concentrate 
on the behavior of the charge as a function of flux ip and 
dot level position £o- 
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FIG. 3: Charge as a function of ip at different values of Ed 
(from top to bottom Ed/F = -0.5, -0.3, -0.1, 0.1, 0.3, 0.5) 
for zero Coulomb interaction (7 = and different A's. Phonon 
frequency Q is much smaller than tunnel resonance width F, 
fi,n/r = 0.05. (a) and (b) Zero Coulomb interaction U — 0. 
A = and the charge Qo is given by Eq. (|16|l at (a). \/Ml — 2 
at (b). The maximum values of the differential charge-to-flux 
sensitivity [Eq. (|18|l ] is near the point 99 = tt with parameter 
£d around zero. For the symmetric barrier and U = the sen- 
sitivity «S(</5, £d) has the singularity at (tt, 0) which disappears 
with any finite asymmetry. Coulomb interaction, or electron- 
phonon interaction, (c) and (d) The same for finite but small 
Coulomb interaction U = T. The "fiat" doublet region ap- 
pears. The width of this region decreases as A increases: from 
(c) to the (d) plot. 



In the absence of phonons (A ~ 0) the dot charge can 
be found from Eqs. (O and (llOp . For the case when the 
singlet constitutes the ground-state 



Qo,: 



TS' 



U/2 



(£0 + C//2)2 + F2 



(16) 



[we should add the electron charge to this result if we 
remember about the subtraction ^^ 1/2 = 1 which ap- 
pears in the dot Hamiltonian ([2])]. For the case of the 
doublet we find: 



gi = o. 



(17) 



[or Qi — e if we restore the constant term which is sub- 
tracted in the Hamiltonian ^] . Everywhere in this arti- 
cle we keep Qi = for symmetry, but the unit charge is 
well defined (no quantum fluctuations) and has the cor- 
rect physical interpretation. Let us start with a normal 
dot with some well-defined integer charge q — 0, e, 2e. 
Then we connect the superconductors through tunnel 
barriers to superconductors (with Cooper pairs). If the 
charge is odd q = e then it does not "feel" the supercon- 
ductors and the charge remains integer. In the case of 



even initial charge g = 0, 2e it couples with Cooper pairs 
in superconductors and creates the singlet state with frac- 
tional and fluctuating charge (with rms value about e). 

For zero Coulomb interaction U = the doublet re- 
gion is absent and for A = the charge is given by Q2 
fom Eq. (fT6|) everywhere, see Fig. [3l[a). Note, that the 
charge Q2 near the values ((p, Ed) = (""jO) has a narrow 
peak and it changes its sign with e^ . For U = \ — this 
peak corresponds to an infinite "charge-to-phase sensi- 
tivity," see Sec. I VII For asymmetric barriers there is 
no point with infinite slope and the maxima of sensi- 
tivity are reached at two locations around ip — tt, see 
Ref. [Ill Note that this peak is broadened by tempera- 
ture, finite superconducting gap A, Coulomb interaction, 
and electron-phonon interaction. The later is shown in 
Fig.Hb). 

Next, if one now considers nonzero Coulomb interac- 
tion [Fig. |3l^c)], then a "flat" doublet region exists: the 
infinitely narrow peak disappears for any U > and in 
the singlet region the charge is still given by Eq. ([TT]). 

In addition, at the boundary of the doublet and the 
singlet region the charge exhibits jumps (for the finite 
superconducting gap A or temperature T > this jump 
is smeared). Therefore the sensitivity is once again sin- 
gular because one abruptly changes the nature of the 
ground-state upon varying ((/9,£o)- 

Further wc study the combination of the charging ef- 
fects in the dot and the electron-phonon coupling. It 
turns out that the size of the doublet region decreases as 
the strength of the electron-phonon interaction A (more 
precisely, the factor X/Ml) increases. This is displayed in 
Figs.[2i;b),[3i;c) andlll^d). Comparing the Figs. [Sljc) and 
md) it can be seen that the overall topology of the plots 
is the same except for the fact that the reduced doublet 
region persists at [/ 7^ 0. 

In the Appendix VK\ we provide a more detailed expla- 
nation of the properties of the function Q((/3, e-o). 



VI. CHARGE-TO-PHASE SENSITIVITY 

One can measure the charge via a capacitively con- 
nected charge detector, e.g., a single-electron transis- 
tor. It does not "feel" the whole Andreev quantum 
dot charge Q in practice, but it feels some renormalized 
charge acOisQ- The geometrical factor ac comes from 
the properties of the measurement gate (capacitance Cm) 
and other capacitively connected parasitic things around 
(capacitance Co): ac = Cin/(Ci„ + Co). The second fac- 
tor as comes from the dynamical feedback of the charge 
detector. In this article we suppose for simplicity that 
ac = as = 1 keeping in mind that the charge in some way 
is suppressed during measurement procedure (we thus 
study the charge sensitivity unperturbed by detector). 

Let us define the charge-to-phase sensitivity at a given 
point as the derivative 
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FIG. 4: (a) The maximal value of charge-to-phase sensitivity 
[see Eq. (|18p ] at (<^, £d) plane as a function of electron-phonon 
coupling strength A at different Coulomb energies U. (b) and 
(c) The values of (p and Ed at which the maximum of the 
sensitivity is attained. At all plots the [/ = (blue), r/3 
(green), 2r/3 (red), F (cyan). The sensitivity can decreases 
or increases with A. At [/ = the sensitivity goes from non- 
physical infinity value at A = 0; this infinity by asymmetry of 
the dot, finite of A, temperature, etc. fifl/F — 0.05. 



This quantity characterizes the charge response to the 
superconducting phases difference, and, hence, to the 
magnetic flux. It can be useful for a flux a detec- 
tor which is based on measuring the charge in the An- 
dreev quantum dotr^ Note that Eq. p^ coincides with 
the current-to-gate voltage sensitivity S = edl/dso 
[={2e^/K)d'^E/de^d{ifl2)]. 

Consider the structure of the sensitivity as a function 
of the parameters [ip, e^,) and its maxima in these param- 
eters iSinax- In this article we concentrate on the sensitiv- 
ity of the singlet region and we omit the sensitivity due 
to the jumps of the charge at the singlet-doublet border. 

For U = \ ~ Q the sensitivity has a "meaningless" 
large value at (<p,£d) — {''^A'^) which corresponds to the 
narrow peak in the charge, see Fig.[31[a). The interaction 
with the vibrating mode cuts this value and the sensitiv- 
ity decreases with A, which is shown in Fig. IDJa) by the 
top line [the maximum moves away from the point (tt, 0) 
and its new position is shown in Figs.SJb) and|H[c)]. 

Given a finite U the sensitivity initially is totally sup- 
pressed by the existence of the non-sensitive doublet re- 
gion, whereas the maximum sensitivity moves to the bor- 
der of the singlet and doublet regions. Increasing A, the 
sensitivity of the singlet region goes down, but the size of 
the doublet region decreases. The competition between 
these two effects gives us new maxima — lower lines in 
Fig. mja). The effect of the decreasing doublet region 
"wins" when the curve goes up (small A's, the maximal 
sensitivity at the singlet-doublet border); when the size 
of the doublet region is small enough the sensitivity has 
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FIG. 5: (a) Subsystem entropy [Eq. (fT9)l ] at £d = and dif- 
ferent v? as a function of A, with hfl/F = 0.05 and U — 0. The 
entropy increases from zero as the parameter X/hQ increases, 
and saturates at value S — log 2. Top inset: dependence of 
the entropy on dot level Ed at superconducting phase differ- 
ence (/3 = (no dependence), tt/2 and n (strong peak around 
Ed ~ 0) for phonon-electron interaction strength X/hQ, = 3. 
(b) Maximum entropy in the {ip,eo) space as a function of 
X/Ml for [/ = (red), U = V (cyan), and U ^ 2T (brown). 
The maximum entropy is reached for e-o — 0, and values of tp 
which depend on U. 



its maximum inside the singlet region and does not de- 
pend on [/ — curves merge and go down. 



VII. ENTROPY 

The entropy provides a measure of the effectiveness of 
the electron-phonon coupling, and in particular to what 
extend this coupling entangles the electron and phonon 
degrees of freedom. The density matrix of the total sys- 
tem is given by p = |$)(^|, where |^) is the ground 
eigenstate of the Hamiltonian H, see Eq. ([5]). The den- 
sity matrices of the electron and phonon subsystems are 
defined as pei = Trph{/5} and Pph = Trei{/5}, respec- 
tively (here Trph and Trei denote traces over electron and 
phonon degrees of freedom). Given a subsystem density 
matrix, the von Neumann entropy is defined as 

S = -Tr{pph log pph} = -Trjpci log pc\}. (19) 

We consider for simplicity a single phonon mode. Let 



us discuss first the absence of Coulomb interaction {U = 
0); then the system can never be in the doublet state, 
and the electronic basis can be restricted to [|0)oi, 12)01]""" 
and the ground-state can be written in full generality 



I*) = flo |0)ci ® bo)ph + a2 |2)oi ® \P2) 



ph; 



(20) 



where |0)oi and |2)oi are the electronic states, |po)ph and 
|p2)ph are normalized phonon states (which can be ex- 
pressed as linear combinations of the basis phonon states 
|n)ph), and |aoP + |Q'2p = 1- The reduced density matrix 
of the electron subsystem is then 



Pel = 



\ao\ 



a-oa2{PQ\P2 



aS«2(p2bo) 
\a2\' 



(21) 



Two extreme cases are simple and notable. First, for 
(Po|P2) = 1, the density matrix corresponds to a pure 
state with zero entropy; accordingly the wave function 
can be factorized as \'^) — (ao|0)ci + a2|2)ci) (81 |po)ph- 
Second, when (po|P2) — 0, then the density matrix 



is diagonal, with S 



|aorlog|ao|^ - |a2rlog|a2| 



which gives the maximum value S = log 2 when |ao| 
|a2p — 1/2. In the general case, the entropy is S* = 
—p+ logp+ — p- logp_, with the eigenvalues p± = 1/2 ± 
[l/4-|aoP|a2p(l-|(pob2)P)]^/^- As entropy is maximal 
{S — log 2) when p+ = p^ = 1/2, and decreases as the 
difference between p+ and p_ increases, we see that to 
have a large entropy one needs to have {po\p2) as small as 
possible, and |aop = 1 — |a2p as close to 1/2 as possible. 

The behavior of the entropy as a function of the pa- 
rameters A (phonon coupling), Sr, (position of the dot 
level) and F — Fcos((/3/2) is shown in Fig.[5ja). 

The main panel of the figure shows the increase in 
the entropy as a function of A, for different values of if, 
and Eo = 0. The fastest increase is obtained for if = tt 
(F = 0). The entanglement of electrons with the other 
degrees of freedom develops easier, if the two electronic 
levels cross or close to each other. The biggest entan- 
glement/entropy ed If — TT then looks natural, since the 
energies E- and E+ [see Eq. ([9])] coincide at the point 
(tt, 0). The difference of E+ — E- increases with increas- 
ing "distance" from point ((/?, £0) to point (tt, 0), and 
correspondingly away from the point (tt, 0) entanglement 
decreases. Also the increasing of the entropy with A can 
be understood from simple analysis of the nature of the 
electron-phonon coupling: the coupling to the electronic 
levels |0)oi and |2)ci "displaces" the phonon field in op- 
posite directions, thus making the phonon states overlap 
|(po 1^2)1 smaller as A increases, which increases entropy. 
Decreasing ip (thus increasing F) gives a smaller entropy. 
This is due to the coupling between the states |0)ci and 
|2)ci when F 7^ 0; the states \po) and \p2) are then com- 
binations of the two displaced states, which makes the 
overlap |(po|P2)| larger and thus decreases entropy. 

The inset of the figure shows how the entropy varies 
when So is changed, for different values of (p: it has a 
peaked behavior, with a width which decreases sharply 
as (fi gets closer to tt (that is, F to 0); for ip = tt, the 



width of the peak is precisely zero: since the electronic 
levels |0)oi and |2)oi are not coupled for </? = tt, any non- 
zero value of Ed rneans that the ground-state is obtained 
with a single electronic state only (|0)ci or |2)oi), and thus 
entanglement with the phonon field, and entropy, is zero. 
Let us now consider the effect of Coulomb interaction 
{U > 0). As has been shown in previous sections, it cre- 
ates a doublet region. There, the entropy is simply zero. 
When a doublet region exists, the maximal entanglement 
between electron and mechanical subsystems is achieved 
at the border of the singlet/doublet regions. Therefore 
the maximum entropy (in variables ((p,eo)) is obtained 
for So = and (f at the edge of the singlet region. The 
entropy as a function of A for non-zero U is plotted in 
Fig. IH^b). The red curve ([/ = 0) is the same as the red 
curve in Fig.lSfa); for U > (cyan and brown curves), the 
entropy maximum goes down because of the existence of 
the doublet region, but again approaches the asymptote 
S = log 2 when the doublet region disappears at large 
A's, see Eq. ^. 



VIII. CURRENT THROUGH THE ANDREEV 
QUANTUM DOT 

The current / = (2e/ft) dE/d{ip/2) is defined by the 
operator 



/=-- |rsin|[dVt-Kh.c.]. 



(22) 



The correspondent matrix elements (in electronic Hilbert 
subspace) 



2e 



Jui^ = -^-i sm — 



V 



1 





10 



(23) 



Similarly to Sec. |Vl] about the sensitivity, we study 

the critical current /c = iaax{I{ip)} dependence on the 

•p 
coupling strength with the vibrational mode A. The later 
behaves in the same way as the sensitivity except for 
special points at (7r,0), see Fig.[S] 

Starting with A = one finds the value of ip which gives 
the maximal value of the current. It can be located at the 
border of the singlet-doublet region (infinitely close from 
the side of the singlet region) and the critical current 
takes the value 



r(D) 



2e 2 
liU 



IP 
4 



{'- + I) 



1/2 



El 

4 



'2 



1/2 



(24) 



Typically, if the critical current is defined hy Iq , its 
value increases with A, e.g., see the red line and top inset 




FIG. 6: Critical current through Andreev quantum dot Ic as 
a function of electron-phonon coupling A. Two effects com- 
pete: electron-phonon coupling dumps the current (see blue 
and green lines, which correspondent to U = and f/ = F) 
and electron-phonon coupling decreases doublet region, where 
/(^) = (increasing of the red line and inset, at bigger 
Coulomb coupling U — 2T). For large \/Ml parameter the 
current totally suppresses. At this plot So = O.IF. 



in Fig. [B] If the critical current's (p is located deep in 
singlet region its value 



r(M) _ ^ J p: 



{r'.2(.„4 



Ed + TT 



(-4)'} 



1/2 



(25) 



decreases with A (e.g., blue line in Fig. [6]). 

The existence of the doublet region transfers the sys- 
tem to the regime of Coulomb blockade; the electron- 
phonon coupling can transfer the system back to the open 
channel regime but simultaneously it partially suppresses 
the current. 



IX. CONCLUSION 

We considered a quantum dot with mechanical de- 
grees of freedom which is coupled to superconducting 
electrodes in a Josephson junction geometry. As such 
a device can be used, in principle, to measure with great 
sensitivity the magnetic flux;^ our main goal was to ad- 
dress the effect of the phonon coupling to this detection 
scheme. The superconducting gap was assumed to be 
larger than all relevant degrees of freedom such as the 
Coulomb energy and the electron-phonon coupling. In 
this so-called "infinite gap limit," retardation effects as- 
sociated with the coupling to the superconducting elec- 
trodes can be neglected, and observables can be com- 
puted using a truncated Hilbert space for the phonons. 

We computed the charge on the quantum dot and de- 
termined its dependence on the phase difference in the 



presence of phonon coupling and Coulomb interaction. 
This allowed to identify regions in parameter space with 
the highest charge to phase sensitivity, which are relevant 
for flux detection. We found that nanomechanical prop- 
erties significantly affect the behavior of the electron sys- 
tem: charge, transport, etc. In the absence of Coulomb 
interaction, the coupling to the vibrational mode reduces 
the charge sensitivity. On the other hand, when Coulomb 
interaction is present, it reduces (eventually completely 
for large coupling) the electrically insensitive doublet re- 
gion due to Coulomb interaction, and in this way it in- 
creases the charge sensitivity. 

Information about the entanglement properties be- 
tween electron and phonon degrees of freedom was ob- 
tained by computing the von Neumann entropy. For a 
fixed phase difference, and in the absence of Coulomb 
energy, the entropy increases with increasing phonon 
coupling, and eventually saturates. When plotted as a 
function of level position, the entropy displays a peak 
when the level position corresponds to the superconduc- 
tor chemical potential. This peak narrows at the phase 
difference approaches tt. When the Coulomb energy is 
switched on, the entropy is zero in the whole doublet 
region. 

Finally, the study of the critical current showed that 
for weak and moderate Coulomb energy, the current is 
typically reduced as the electron-phonon coupling is in- 
creased. For a larger Coulomb coupling which exceeds 
the dot line width, the critical current is much reduced 
at small electron-phonon coupling but it acquires a max- 
imum for larger coupling strength and eventually merges 
with the curves corresponding to weak Coulomb interac- 
tion. 
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FIG. 7: Charge as a function of {eu, f) for zero U — [(a)-(c)] 
and nonzero U — 2r [(d)- (f)] Coulomb interaction and differ- 
ent A's. Single phonon mode with frequency Q much smaller 
then tunnel resonance width F, ftfi/F — 0.05. (a) A = 0, 
the charge Qo is given by Eq. (|16|) . The maximal values of 
differential charge-to-flux sensitivity [Eq. (jlSp ] lies near point 
{if — TTjErt =0). For a symmetric barrier and [7 = sensitiv- 
ity S{ip,eu) has the special point at (tt, 0) but it disappears 
with any finite asymmetry. Coulomb interaction, or electron- 
phonon interaction, (b) With increasing vibration coupling 
\/Ml = 1.5 the maximal charge-to-flux sensitivity decreases, 
(c) \/Ml — 3. (d) Nonzero Coulomb interaction U = 2r, no 
electron-phonon coupling A = the largest doublet region is 
described by Eq. ([TT]). (e) Increased \/Ml = 1.5 decreases the 
doublet region (f) X/Ml — 3. For a finite temperature T the 
border of the single-doublet regions is smeared, with a width 
of fcBTin (d)-(f). 



Appendix A: Detailed explanation of the charge 
behavior 



For zero Coulomb interaction U = the doublet region 
is absent and for A = the charge is given by Eq. p^ 
everywhere. For U = X = the derivative dQ/df has 
a jump at the location {if,e-o) = (tt, 0) in the 2D plane 
{ip^Eo), see Fig.ITJa); near this point the charge-to-phase 
sensitivity tends to infinity. Note that for asymmetric 
barriers this point with an infinite slope does not exist 
and maxima of sensitivity are reached at four locations 
around (tt, 0), see Ref . Illl. When the electron-phonon in- 
teraction is switched on this special point disappears and 
the maximal sensitivity is therefore suppressed by the 
electron-phonon interaction. This is displayed with in- 
creasing electron-phonon coupling in Figs.jT^b) and[7l^c). 

Next, if one now considers nonzero Coulomb interac- 
tion [Fig. [Tljd)], then a "flat" doublet region exists: the 



special point with infinite sensitivity disappears for any 
U > Q and the sensitivity in the neighborhood of this 
point is totally suppressed by the existence of the dou- 
blet region. In the singlet region the charge is given by 
Eq. (HID. 

In addition, we observe that at the boundary between 
the doublet and the singlet region, the charge exibits 
jumps and therefore the sensitivity is once again singu- 
lar (for finite temperatures T > this jump is smeared) 
because one abruptly changes the nature of ground-state 
upon varying {tp,£v)- 

We next study the combination of electron interaction 
on the dot with the electron-phonon coupling. By in- 
creasing (from zero) the strength of the electron-phonon 
interaction A, the size of the doublet region decreases. 
This is displayed in Figs. [7{e) and [71[f ) where the same 
electron-phonon coupling parameters are chosen as in 



Figs.[7Ub) andlTjc). The evolution of the doublet region 
is plotted in Fig.[2]for more values of the electron-phonon 
coupling constant. The largest loop corresponds to the 
smallest (zero) A. Upon switching A, the reduction of 
this region is barely noticable, it acts mostly on the level 
position range as it still approaches the phase values 
and 27r. Beyond A = 0.4, the reduction is effective in 
both the {(p^e-o) direction. 



There is a competition between charge repulsion ef- 
fects on the dot and the presence of the electron-phonon 
coupling, which can be understood to be playing the role 
of an effective attractive, retarded, interaction. This ex- 
plains the reduction of the doublet region. When compar- 
ing Figs.[71[c) and[71[f) we note that the overall topology 
of the plots is the same except for the fact that a reduced 
doublet region persists at C/ 7^ 0. 
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